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Particle decay and 21 cm absorption from first minihaloes 
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ABSTRACT 

We consider the influence of decaying dark matter (DM) particles on the character- 
istics of 2f cm absorption in spectra of distant radio- loud sources - "21 cm forest" - 
from minihaloes with masses M = 10 5 - 10 7 Mq virialized at z mr = 10. We use ID 
self-consistent hydrodynamic description to study evolution of minihaloes, and follow 
up their absorption characteristics from turnaround to virialization. We find that in 
the presence of decaying dark matter both thermal and dynamical evolution of mini- 
haloes demonstrate significant deviation from those in the model without dark matter 
decay (standard recombination). We show that optical depth in 21 cm line is strongly 
suppressed in the presence of decaying particles: for M = 10 5 — 10 6 Mq decaying dark 
matter with the energy rate deposited in baryonic gas £l = 0.59 x 10~ 25 s _1 - the cur- 
rent upper limit of the energy deposit - decreases the optical depth and the equivalent 
width by an order of magnitude compared to the standard recombination. Thus addi- 
tional ionization and heating from decaying DM particles almost "erases" absorption 
features from minihaloes with M = 10 5 — 10 6 Af Q for £ £ 0.3£l, which consequently 
considerably decreases the number of strong absorptions: for example, the number of 
absorptions with the equivalent width W° hs Z, 0.3 kHz at z ~ 10 decreases more than 
2.5 times for = 0.3 and £ 4.5 times for = 1. We argue that "21 cm forest" 
absorptions might be a powerful probe of the presence of decaying dark matter in the 
early Universe. 
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1 INTRODUCTION 

A way to study pre-reionization history of the inter- 
galactic medium lies in observation of the redshifted 21 cm 
line produced in the hyp e rfine transition of neutral hydro- 
gen (iHogan fc Reedfl~979l : IScott et al. Ill99ll : iMadau et aT\ 
1 19971 ). There are at least three techniques to use the 21 cm 
line for studying this epoch: observations of emission or ab- 
sorption of the neutral intergalactic medium (IGM) against 
the cos mic microwave background (CM B) - the 21 cm global 
signal ll Shaver et al. 1999: Sethi 2 0051 ). statistical studies of 
the angular distributio n of the intensity i n the 21 cm line - 
the 21 cm fluctuations (|Tozzi et al. Il2000l : llliev et al. |2002h . 
and measurements of absorption along the lin e of sight to a 
distant radio-loud quasar - the "21 cm forest" (|Carilli et al. I 
l2002l : iFurlanetto fc Loebll2002T ) . The first two methods pro- 
vide mapping of the 3D space-redshift features of the IGM, 
whereas the latter works only along line of sight, i.e. pro- 



vides a ID probe. A physical distance corresponded to the 
angular resolution of existing and future radiotelescopes, 
such as e.g. lofarQ, is about one megaparsec. There- 
fore only huge structures, such as for example, HII regions 
formed by stellar clusters and quasars, the large scale cos- 
mic web structure, are possible to detect. Instead, by mea- 
suring the absorption features along redshift we can resolve 
small scale structu re, e.g. individual minihaloes, galaxies , 
stellar HII regions llKumar et fflT1ll995l: iBagla et al. I | l997| : 



Furlanetto fc Loebll2002l : lFurlanettdl2006t lYue et al. II2OOS 
Mack fc Wvithell201ll ). 

Ionizing photons produced by first stellar objects 
change gas thermal state in minihaloes and their vicin- 
ity, and as a consequence, affect the 21 cm absorption 
features (|Yue et al. 1120091 ; IXu et al. 1120111 ). Similar effect 
can be expected from other possible sources of ioniz- 
ing photons, such as decaying dark matter (DM) parti- 
cles (|Sciamal 1 19821 : IScott et al\ Il99ll : iBharadwai fc Sethil 
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Chen fc Kamionkowskill2004lHansen fc Haimanll2004l; 
Kasuva fc Kawasakil |2004|; iKasuva et al. I 20041 ; IPierpaolil 
2004 iBelikov fc Hooper! I2009T), annihilation of dark mat- 



ter dChuzhovl 120081; iMvers fc Nusseil 120081 : lYuan . et al. 1 
l2010l ; ICumberbatch et al. II2OIO1 ) and ultra-high energy cos- 
mic rays produced in decay of superh eavy dark matter 
parti cl es with masses M x £ 1 12 GeV dBerezinskv et ai] 



19971; iBirkel fc Sarkarl 1 19981; 



Doroshkevich fc Naselskvl 12002] : 



Kuzmin fc Rubakov 



Doroshkevich et 



1908: 
l2003h . 



Previous studies have shown that ionizing background pro- 



evolution of the intergalactic medium jDodelson & Jubat 


1994 


; Biermann & Kuscnko 


20061; Shchekinov & Vasiliev 


2004 


;l Vasiliev & Shchekinov! 


20061; iMapelli et al. 20061") and 



influe nce through it the 21 cm global signal and fluctua- 
tions dFurlanetto et aFlbOOrj; IShchekinov fc Vasilieyl p007 



Chuzhovl 120081; IMvers fc Nusseil 120081; lYuan et al I boiC 



Cumberbatch aZ. 1 12010| : iNataraian A. fc Schwarzj 12009? ) 
Recent studies stron gly constrain properties of such ioniz- 
ing radiation sources dZhang et al. ll2007l;|Cirelli et al. 120091: 



DeLope Amigo et al 



2009 



ng et at. II^UUYNUirelli et al. Iz UUa ; 
9UPeter e£ oTtoloHZhang ei al. \ 



20101 : iGalli et al. Il201ll ; iHiitsi et aTll201ll ), though do not 



fully exclude them. Therefore, in this paper we focus on how 
the 21 cm absorption - "21 cm forest" - depends on the de- 
caying dark matter parameters. 

The 21 cm forest is supposed to be a signal from nu- 
merous minihaloes along line of sight being at different 
evolutionary states from the turnaround to the virializa- 
tion, with the absorption features depending on the evo - 
lutionary state |Meiksinll20"Tll ; IVasiliev fc Shchekinovll2012l ). 
In several previous studies minihaloes were considered as 
static objects with fixed dark matter and ba ryonic pro- 
files (|Furlanetto fc Loebl |2002| ; lFurlanettdl2006|'). More re- 



cently accretion has been added (|Xu et al. I 1201 j ). and 
further, the evolutionary effects during th e format i on of 
minihaloes have been taken into account (|Meiksinl l201ll ; 
IVasiliev fc Shchekinovll2012l ). 

Another aspect important when characteristics of the 
21 cm forest are concerned relates to spatial distribution of 
matter in min ihaloes. The dark ma tter profile is mostly as- 
sumed cuspy (|Navarro et al. |[l997f ). However, according to 
the observati ons of local dw arf galaxies the dark matter pro- 
file is flatter |Burkcrt 1995|). Recently in several theoretical 
studies it has been pointed out that th e first protogalaxies 
may also have a flat dark matter profile ( Mashch enko et al. I 
120061 ; iTonini et al. 1 120061 ; iMikheeva et al. 112007ft . It seems 
obvious that changes in dark matter profiles result in cor- 
responding changes of the radial distribution of baryonic 
component inside minihaloes, and consequently, in the 21 
cm absorption properties. 

In this paper we study the effects of decaying 
dark matter particles on absorption features of non- 
static (evolving) minihaloes with a flat dark matter 
profile. We assume a ACDM cosmology with the pa- 
rameters (Q , f2 A , n m , tl b , h) = (1.0,0.76,0.24,0.041,0.73) 
|Spergel et al. 1120071 1. 



2 EVOLUTION OF MINIHALO 

To model the evolution of minihaloes we use a one- 
dimensional Lagrangian code. For dark matter we follow the 



description given by iRipamontil |2007l ). Namely, the dark 
matter, Mdm = ^DAiMhaio/^M, is assumed to be enclosed 
within a certain truncation radius Rt r , inside which the 
dark matter profile is a truncated isothermal sphere with 
a flat core of radius R cor£ . The parameter 77 = R CO re/ Rvir 
is taken 0.1 for all simulations. Such a description is used 
to mimic the evolut ion of a simple top-hat fluctuation (e.g. 
|Padmanabhanlll993l '). In order to describe time evolution of 
the truncation radius we make use of the approximation for- 
m ula for a top-hat density perturbation in the form proposed 
bv lTegmark et al. I (| 19971 ). 

For simulations of baryonic component we use a 
ID Lagrangian s c heme similar to that described by 
iThoul fc Weinberg! (|l995l ). As a standard resolution we used 
1000 zones over the computational domain, and found a rea- 
sonable convergence. 

Chemical and ionization composition include a stan- 
dard set of species: H, H + , H~, He, He + , He ++ , H2, H^", 
D, D + , D~, HD, HD + and e. The corresponding reac- 
tion r ates are taken from (|Galli fc Palla 1998: St ancil et al. I 
1998!). Energy equation includes radiative losses typical in 
primordial plasma: Compton cooling, recombination and 
brem sstrahl ung radiation, collisi onal excitati on of H I (ICenl 
19921). H? (IGalli fc Palla! Il998h and HD jFloweri I2OO0I ; 



Lipovka et al. 1120051 ) 



We start our one-dimensional simulations at redshift 
2 = 100. The initial parameters: gas temperature, chemical 
composition of gas and other quantities - we have taken 
from simple one-zone calculations began at z = 1000 with 
typical values at the end of recombination: T gas = Tcmb, 
x[U] = 0.9328, x[U + ] = 0.0672, x[D] = 2.3 x 10~ 5 ,a;[p + T = 
1.68 x 10 -6 (see references and details in IRipamontil I20071 
Table 2). 



3 IONIZATION AND THERMAL HISTORY 

The contribution from decaying dark matter is de- 
scribed by the corresponding ionization and heating terms 
into equations for the ionization and thermal evolu- 
tion. The ionization rate due to presence of additional 
ionization sources from deca ying particles is written as 
(|Chen fc Kamionkowskll2004h 



hv c 



(1) 



where ~Xi is the energy frac tion deposited into ionization 
(IShull fc van Steenbergj[T985l ). m v is the proton mass, f x = 
Qx(z)/£lb(z), £lb(z), the baryon density parameter, ilx(z), 
the fractional abundance of decaying particles, Fx is the 
decay rate, hu c , the energy of Ly-c photons. In general, the 
ionization history depends on the energy rate deposited in 
baryonic gas, £ = Xifx^x (|Chen fc Kamionkowskil [20041 ). 

The heating rate produced by ionizing photons from 
decaying dark matter particle s can be written in the form 
l|Chen fc Kamionkowsklll2004h 



K = XhrripC f x T) 



(2) 



where ~Xh is the energy fra ction depositing into heating 
(|Shull fc van Steenberg|[l985l ). By order of magnitude \i ~ 
Xh ~ 1/3 for the conditions we are interested in. 
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Using the CMB datasets lzhang et a~L~\ (|2007h have con- 
strained the energy rate of the radiatively decaying dark 
matter as the upper limit £l 1.7 x 10 -25 s" 1 . Extend- 
ing analysis with data of Type la supernova, Lya forest, 
large scale structure and weak l ensing observations have 
lead iDeLope Amigo et al. I (|2009l ) to stronger constraints: 



6 



0.59 x 10" 



It is worth noting that all datasets 



favour long-living decayi ng dark matter particles wi th the 
lifetime F" 1 > 100 Gyr (|DeLope Amigo et al. 1 120091 '). Fur- 
ther improvement is expected from the Planck satellite. 
In this paper we consider the decaying dark matter with 
£ ^ £l = 0.59 x 10~ 25 s _1 and compare results with the 
standard recombination scenario, i.e. £ = 0. 



4 SPIN TEMPERATURE 

The two processes: atomic collisions and scattering of ultra- 
violet (UV) photons, couple the HI sp i n temperature and 
the gas kinetic temperature (|FieIdlll95a : IWouthuvsenlll952T ) 



T s = 



Tcmb + y a Tk + ycTk 



(3) 



1 + Va + Vc 

here Tcmb is the CMB temperature, y c , y a are the functions 
determined by the collisional excitations and the intensity 
of the UV resonant photons 



PwT* 
AwT k 



Vc = 



C 10 T. 
AioT k 



(4) 



T, = 0.0682 K is the hyperfme energy splitting, Aio = 
2.87 x 10~ 15 s _1 is the spontaneous emission rate of the 
hyperfine transition, Cio = kionn + "/en e is the collisional 
de-excitation rate by hydrogen atoms and electrons, the pho- 
ton induced de-excitation rate is neg ligible , for fcio we use 



the approximation bv lKuhlen et al. (|2006h . for j e we take 
the approximation from LisztJ 1 200 ll ). Pio is the indirect de- 
excitatio n rate, whi ch is related to the total Lya scattering 
rate P a (|Fieldlll958r ) 



Pio = 4P a /27, 



(5) 



where P a = J cn v a u dv, n„ is the number density of photons 
per unit fr equency range, a(v) is the cross-section for Lya 
scattering (|Madau et al. II 19971 ). 

Additional ionizations produced by decaying dark mat- 
ter deposit into Lya photon budget due to recombination. 
In our calculations we do not take into account the UV 
pumping, firstly because such deposit is minor in compar- 
ison with the atomic collisions for the decaying dark mat- 
te r parameters and redshifts c onsidered here (see Figure 2 
in lShchekinov fc Vasilievll2007l ). Secondly, contribution from 
Lya background produced by star-forming minihaloes is 
small due to the rarity of first stellar objects at redshifts 
2 = 10— 15. Moreover, it can be readily shown that in order 
to compete contribution from atomic and electron collisions 
to HI spin temperature UV energy flux has to be unaccept- 
ably high: F Lya > 0.1 erg s _1 cm -2 for gas with T ~ 200 K 
and n ~ 1 cm -3 at z ~ 10. 



5 RESULTS 

Stellar and quasi-stellar sources of ionizing radiation be- 
gin to form at redshifts z < 20. They heat the surround- 



ing gas through photoionization, which then emits in 21 
cm with a patchy, spot-like distribution on the sky. Instead, 
the decaying particles illuminate and heat the IGM homo- 
geneously by X-ray photons normally produced in decays 
jChen fc Kamionkowskill2004l ). Such photons have mean free 
path length much longer than the typical distance between 
minihaloes, d ~ 30(1 + z/20) kpc. Therefore in minihaloes 
the ionization and heating rates due to decays of dark mat- 
ter particles can be taken equal to the rates in the IGM 
background. 



5.1 Dynamics of minihaloes 

Heating from decaying dark matter weakens the accretion 
rate of baryons on to the minihaloes. Therefore, in compari- 
son with the evolution of minihaloes in the standard scenario 
one can expect smaller baryon mass and higher temperature 
within the virial radius of a minihalo. The influence of de- 
caying particles depends on the minihalo mass and the rate 
of the decay energy deposited in gas, £. 

Here we consider evolution of minihaloes within the 
mass range M — 10 5 — 1O 7 M0 virialized at z v = 10. Our 
choice of the redshift of virialization is caused by the fact 
that at lower redshifts the effects of reionization begin to 
play dominant role and thus contaminated contribution from 
decaying dark matter. On the other hand, at higher red- 
shifts the number of bright radio objects which can serve as 
background sources where minihaloes can imprint becomes 
low. In the mass range we focus on the three character- 
istic masses: M = 10 s , 10 6 and 10 7 M Q for the follow- 
ing reasons. It has been found that inside minihaloes with 
M ~ 2 x 10 6 Mq the star formatio n can potentially occur 
at z ~ 10 (see also iRipamontillboOTT ) ■ in the sense that such 
minihaloes collapse at z V i r — 10, i.e. the gas density in the 
most inner shell reaches around 10 s cm~ 3 and grows up fur- 
ther. Minihaloes with lower masses cannot cool efficiently 
and form stars at z ~ 10. Therefore the first mass considered 
here M = 1O 5 M0 represents in general the low-mass end of 
non-starforming minihaloes. The second value M = 10 6 Mq 
cannot also collapse at z ~ 10, but it is close to the critical 
mass M ~ 2 x 1O 6 M0, and represents qualitatively correct 
characteristics of 21 cm absorptions in the range of masses 
of marginally star-forming minihaloes. The third value is 
M = 1O 7 M0: such minihaloes collapse (and therefore can 
form stars) even before the virialization redshift, namely, at 
2 ~ 11. Thus, we study evolution and absorption charac- 
teristics of the minihaloes, which from one side can remain 
dark after virialization (M = 10 s — 1O 6 M0) and from the 
other form first stars (M = 1O 7 M0). 



5.1.1 Evolution of minihaloes with M — 1O J M0 and 
10 6 M o . 

Figures [1121 show radial profiles of number density (upper), 
temperature (middle) and velocity (lower panels) in mini- 
haloes M = 10 5 , 1O 6 M0 virialized at z V i r = 10, for the three 
redshifts z = 15.5, 12 and 10 (from left to right panels, cor- 
respondingly) . 

For M = 10 5 M Q the virial temperature about 500 K is 
close to the one in the IGM heated by decaying DM parti- 
cles with ^, 0.1. At the turnaround state (left panels) 
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Figure 1. The radial density (upper), temperature (middle) and velocity (lower panels) profiles of a halo M = 10 Mq virialized at 
Zvir = 10 are shown at redshifts z = 15.5, 12 and 10 (from left to right panels, correspondingly) in the standard recombination model 
(dots) and in the presence of decaying dark matter with = 0.1,0.3,1 (solid, dash and dot-dashed lines, correspondingly), where 

f L = 0.59 X 10 -25 b _1 . Note that in the presence of decaying dark matter gas almost ceases collapse into the minihalo potential well. 




Figure 2. The radial density (upper), temperature (middle) and velocity (lower panels) profiles of a halo M = 10®Mq virialized at 
Zvir = 10 are shown at redshifts z = 15.5, 12 and 10 (from left to right panels, correspondingly) in the standard recombination model 
(dots) and in the presence of decaying dark matter with = 0.1,0.3,1 (solid, dash and dot-dashed lines, correspondingly), where 

£l = 0.59 X 10 -25 s — 1 . Note the absence of the accretion (virial) shock for 0.3 (right lower panel), whereas for the lower 

the profiles tend to those in the standard recombination model. 
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Figure 3. The radial density (upper), temperature (middle) and velocity (lower panels) profiles of a halo M = 10 7 Mq virialized at 
Zvir = 10 are shown at redshifts z = 15.5, 12 and 11 (from left to right panels, correspondingly) in the standard recombination model 
(dots) and in the presence of decaying dark matter with = 0.1,0.3,1 (solid, dash and dot-dashed lines, correspondingly), where 

£l = 0.59 X 10 — 25 s" 1 . Note that the velocity profiles demonstrate a stable gas collapse in the potential well of the dark halo. The 
density and velocity profiles for all considered arc close to those in the standard scenario, though the temperature profiles 

demonstrate appreciable difference. 



temperature profiles differ significantly from the standard 
one: namely, they show the internal and the external tem- 
peratures practically equal due to additional heating from 
decaying DM particles. Such a high background temperature 
in the presence of decaying particles can cause stopping gas 
collapse in the minihalo potential well during further evolu- 
tion. 

At z = 12 (middle panels of Figure [TJ the standard 
model shows obvious gas accretion resulting in an increase 
of density and temperature. At the same time in models with 
£ Si 0.3£l gas velocities are non-negative at any radius, so 
that the baryon content of minihalo does not grow. The den- 
sity within the plateau r ^ 100 pc for £ = £l is an order of 
magnitude lower than in the standard recombination - this 
is clearly seen on the lower panel of Figure [T] At the virial- 
ization this difference in density increases: for low ionization 
rates £ £ 0.3£l the density grows several times compared to 
that at z — 12, while for higher £ = £l it stays nearly the 
same. As a result for high ionization rates the deviation from 
the standard model reaches two orders of magnitude. Practi- 
cally in the whole range of ionization rates temperature pro- 
files within r 100 pc do not deviate significantly from each 
other, showing a dominance of the decaying DM in heating. 
The most appreciable difference between the standard re- 
combination and the models with decaying dark matter is 
the absence of a significant accretion shock wave around 
r ~ 100 pc. In the standard case such a wave is obviously 
seen on the right lower panel of Figure [1] Thus the presence 



of decaying dark matter prevents collecting baryonic mass 
in the minihaloes with M ~ 1O J M0. 

For £ XL 0.3£z, the evolution of minihaloes with M = 
10 6 M Q is similar to that with M = 10 6 M© (the dash and 
dash-dotted lines in Figure [2]). Whereas for lower rate £ evo- 
lution of M — lO 6 A/0 looks closer to the standard recombi- 
nation model. Note, in particular, the presence of the virial 
shock for £ = 0.1£l (solid lines at the lower row of panels in 
Figure [2]) similar to that in the standard case (dotted line). 

5.1.2 Evolution of minihaloes with M = 1O 7 M0. 

For M — 1O 7 M0 the density profile^ for all £ considered 
here are close to those in the standard recombination sce- 
nario (Figure [3}. Instead, temperature and velocity profiles 
demonstrate considerable differences. These differences can 
be explained by the accretion of warmer intergalactic gas. At 
the turnaround temperature inside the halo is almost equal 
to the background value (middle left panel), which is higher 
than the standard background temperature. At this stage 
gas has already begun collapsing into the potential well of 
the dark halo, providing the velocity profiles almost coin- 
cided for all £ £l (lower left panel). In further accretion 
of warm intergalactic gas the differences between models 

2 As for minihaloes with M = 10' Mq the collapse (i.e. the time 
when gas density in the most inner shell reaches 10 s cm -3 and 
grows further) occurs earlier than formal virialization, wc present 
the results for not later than the collapse redshift, i.e. for z = 11. 
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Figure 4. Optical depth at the impact parameters a = 0.1, 0.3, 1, 1.5, 3r/r„; r (solid, dash, dot, short dash and dot-dash lines, 
correspondingly) for a halo M = 10 s Mq virialized at z v i r = 10, at redshifts z = 15.5,12 and 10 (from top to bottom panels, corre- 
spondingly) in the standard recombination model and in the presence of decaying dark matter with = 0.1, 0.3, 1 (from left to right, 
correspondingly), where = 0.59 X 10 — 25 s" 1 . 



start revealing at z = 12. Temperature inside the minihalo 
(r ^ 500 pc) becomes considerably higher for the models 
with decaying particles. The absolute value of the velocity 
decreases with an increase of the rate f, i.e. baryons col- 
lapse slower. At the collapse redshift, z ~ 11, differences in 
the temperature and the velocity profiles for different rates 
£ remain significant in the whole radial range. In all mod- 
els there the virial shock wave does not appear, the velocity 
profiles demonstrate stable collapse of gas in the potential 
well of the dark halo. One can note that the profiles for the 
models with decaying particles are close each other and to 
the profile for the standard reionization within r 200 pc 
which corresponds to r £ 0.5r v i r . 

The evolutionary differences for minihaloes in the mass 
range M = 10 s — 1O 7 M0 in the decaying dark matter models 
are expected to become seen in the HI 21 cm absorption 
characteristics. In the following sections we consider how the 
differences manifest in optical depth and equivalent width 
of the 21 cm absorption. 



5.2 Optical depth 

The optical depth at a frequency v along a line of sight is 



3h p c 3 A 10 



d:v 



riHi(r) 



exp 



y^b 2 (r)T 3 (r) 
v(v) - vi(r)] 2 



(6) 



b 2 {r) 



where r 2 = (ar„, r ) 2 + x 2 , a = r/r V i r is the impact pa- 
rameter (a line of sight crosses the minihalo at a distance 
ar V i r from the centre of the minihalo), v{v) = c{v — uo)/vo, 
vi(r) is the infall velocity projected on the line of sight, 
b 2 = 2kTk(r)/m p is the Doppler parameter of gas. 

Figures |4][6] show the optical depth at several impact 
parameters, a, for the haloes with M = 10 5 , 10 6 , 1O 7 M 
at redshifts z = 15.5, 12 and 10 (from top to bottom, cor- 
respondingly jf| in the standard recombination model and in 
the presence of decaying dark matter with = 0.1, 0.3, 1. 
The analysis of the 21 cm absorption features of mini- 
hal oes in the standa r d recombination scenario can be found 
in JXu et al I 1201 ll ; iMeiksinl l201ll ; IVasiliev fc Shchekinov! 
l2012i ). 



5.2.1 Line profiles in on-star-forming minihaloes: 
M = 1O 5 M and 1O 6 M . 

Minihaloes with M = 10 s Mq show a small peak in optical 
depth at v — 1/0 — 40 kHz (FigureU]). This peak emerges from 
dynamics of the halo: it forms in a narrow external layer of 
the halo where a non-zero diverging velocity overlaps with a 
fast decreasing temperature profile - b 2 (r) in the denomina- 
tor of the exponent (|6]). The peak weakens in models with 
decaying DM. Note very weak or lack of dependence of the 
optical depth at v — vo > 10 kHz on the impact parameter, 
reflecting the fact that the velocity profiles are flat within 



A For M = 10 'Mq the lowest redshift is z 
for details. 



11, see Section 5.1 
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Figure 5. The same as in Figure[4] but for a halo with M = 10 6 Mq: a = 0.1, 0.3, 1, 1.5, 3r/r v i r (solid, dash, dot, short dash and 
dot-dash lines, correspondingly), at redshifts z = 15.5, 12 and 10 (from top to bottom). 
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Figure 7. The observed equivalent widths for minihaloes with mass M = 10 s , 10 6 , 10 7 M Q (from top to bottom) at redshifts z = 15.5, 12 
and z = z v i r = 10 (dot-dash, dash and solid lines, correspondingly) in the standard recombination model and in the presence of decaying 
dark matter with = 0.1, 0.3, 1 (from left to right, correspondingly), where £l = 0.59 X 10 — 25 s _1 . Note that the y-axes are different 

between panels. 



haloes and do not practically depend on the ionization pa- 
rameter {TJ. I n the standard scenario optical depth drops 
sharply at higher frequencies, v — Vo > 50 kHz, from one 
side because temperature at the boundary falls down to very 
low background temperature, and from the other, because 
velocity at the boundary joins the Hubble flow (such that 
the numerator in the exponent of © increases). In the pres- 
ence of decaying dark matter the peak at v — Vq ~ 40 kHz 
is suppressed due to the increase of the background temper- 
ature (leading to an increase of denominators both in the 
exponent and in the integrand of (©)■ For higher masses 
similar peaks are found at higher frequencies, though they 
contribute little to absorption. 

An obvious result is that the presence of decaying dark 
matter particles heavily suppresses optical depth: for M = 
10 s Mq in the standard recombination scenario optical depth 
at a = 0.1 reaches ~0, and falls below ~0.06 for £ > 0.1£l- 
For M — lO 6 A/0 optical depth shows similar dependence, 
though frequency profiles at a > 0.3 seem in principle dis- 
tinguishable at lower ionization rates £ 0.3£l (Figure [SJ. 
Thus, the decaying dark matter practically erases 21 cm ab- 
sorption features from minihaloes with M = 10 5 — 1O 6 M0 
and the ionization rate higher than £ £ 0.3£l- 

Horn-like profiles in M = 10 6 Mq haloes at z = 12 (mid- 
dle left panel of Figure [5]) are due to continuing accretion of 
atomic hy drogen at t h ese st ages - such profiles are first de- 
scribed by IXu et all |201ll ). They disappear either during 
further evolution when the halo settles down to a virial state 
(lower left panel), or in the models with an increase of the 
decaying dark matter where heating prevents contraction 
(middle row of panels). 



5.2.2 Line profiles in star-forming minihaloes: 
M = 1O 7 M . 

In the standard recombination horn-like profiles due to 
strong accretion i n massive haloes around its virialization 
|Xu et al. I l201ll ; IVasiliev fc Shchekinov! |2012| ) are clearly 
seen at a ^ 0.3. More precisely, they origin due to lower spin 
temperature in the region with the maximum infall velocity 
as seen from the integrand in (JUJ). Higher infall velocity for a 
minihalo M — 1O 7 M0 shifts maximum of the optical depth 
towards higher frequencies (cf. the position of the optical 
depth maximum at z — 12 and z — 10 for e.g. a = 0.1, see 
Figure [3| . 

In general, in the presence of decaying dark matter the 
accretion rate (i.e. the infall velocity) weakens (see Figure^, 
because of warming of the accreting gas. As a consequence 
the optical depth maximum shifts to lower frequencies and a 
horn- like dependence of the optical depth almost disappears 
for££0.1£L. 



5.3 Equivalent width 

The suppression of the 21 cm line optical depth by decay- 
ing DM particles manifests in decrease of the observed line 
equivalent width determined as W° bB = W v /(1 + z), where 
the intrinsic equivalent width is 



')dv 



(1 



')dv 



{-) 



where tigm is the optical depth of the background neutral 
IGM. Figure [7] demonstrates the observed equivalent width 
for minihaloes with mass M = 10 5 , 10 6 , 1O 7 M at redshifts 
z — 15.5, 12 and 10, for the standard recombination and in 
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the presence of decaying dark matter with = 0.1, 0.3, 1. 
As expected the equivalent width decreases significantly 
with an increase of £: the equivalent width at the impact 
parameter r/r„i r ^ 0.5 for M = 10 J Af© at z = 10 (solid 
lines) decreases from ~ 0.2 KHz for the standard recombi- 
nation to a negligible value 0.01 kHz for = 1. At 
2 = 15.5, 12 (dash and dot-dash lines) the equivalent width 
shows not such a catastrophic drop, however it hardly can 
be resolved on future telescopes a nd deposit into the total 
signal in broad-band observations |Xu et al. 1120111 ). 

Lines with equivalent widths W° bs £ 0.2 kHz can be 
formed by haloes with M > 1O 6 M at z = 10 (solid lines) 
for ~ 0.3 up to a — r/r„i r 1. In the higher redshift 
range, z = 15.5, 12, such strong lines can be associated only 
with M ^ 10 6 Mq in the standard recombination model. 
Higher equivalent widths ~ 0.5 — 0.9 kHz can be attributed 
to low impact factors a ~ 0.3 in the standard recombina- 
tion. However, more massive haloes, M = 10 7 M©, can form 
strong absorptions with the width £ 0.5 kHz even for high 
ionization rate £ ~ £&, though only in a narrow impact fac- 
tor range a £j 0.1. Thus, the decaying dark matter leads to 
a considerable decrease of the number of strong absorption 
lines. This suppression of the 21 cm absorption is similar 
to what stems from the action of the X -ray background 
jFurlanetto fc Loebll20o3 ; IXu et al. Ilioill ). For £/£ L > 0.3 
strong lines are formed only in high-mass (M ~ 1O 7 M ) 
haloes nearly along the diameter, but such events are rarely 
expected. 

The sensitivity of future telescopes is more than one 
order of magnitude lower than the flux limit needed to 
separate spectral lines from individual haloes, i.e. to ob - 
serve with spectral resolut ion of 1 kHz (|Xu et al. 1 1201 if ) . 
Therefore, IXu et al. I |201ll ) have proposed broad-band ob- 
servations with lower resolution. In this case one measures 
averaged signal from different haloes crossed by a line of 
sight. Accordingly, we introduce a radially averaged ob- 
served equivalent' width: (W° bs ) = J^ mr W° bs (r)rdr /r 2 vlr . 
Figure [8] presents the dependence of this value for mini- 
haloes M = 10 s , 1O 6 M , 10 7 M o on The increase 
of (W° bs ) for M = 10 5 M Q at z = 10 and \/£ L < 0.1 is ex- 
plained by a widening of W„ bs in the range of impact factor 
up to r/r V i r 2 (see upper panels in Figure [7J. For other 
considered masses and redshifts the radially averaged W° bs 
decreases when grows by a factor of 1.2-2, which appar- 
ently is sufficient to identify effects from additional ioniza- 
tion sources, and in particular, from decaying DM particles. 

In previous studies theoretical spectrum of the 21 cm 
forest has been simulated within the assumption of steady- 
state mini haloes with fixed profiles c orresponding to the viri- 
alization jFurlanetto fc Loebl [20021 ; IXu et al. Il201lf ). When 
dynamics is taken into account, minihaloes lying on a given 
line of sight are, in general, at different evolutionary stages, 
and have different density, velocity and temperature pro- 
files. Moreover, Press-Schechter formalism cannot anymore 
be used for description of their mass function as soon as 
they are lying far from the virialization. It makes problem- 
atic to describe correctly the number density of minihaloes 
at each evolutionary stage. Possible formation of stars in 
massive evolved minihaloes brings additional complication. 
Overall, the analysis of theoretical spectrum within a sta- 
tistical simulation becomes exceedingly cumbersome. More 
relevant picture can be obtained from high-resolution cosmo- 



A 



N 




Figure 8. The dependence of the radially averaged ob- 
served equivalent width for minihaloes with mass M = 
10 5 , 1O 6 M0, 1O 7 M0 (from top to bottom panels, correspond- 
ingly) at redshifts z = 15.5, 12 and 10 (dot-dash, dash and solid 
lines, correspondingly) on the normalized rate of the decay energy 
deposited in baryonic gas where = 0.59 X 10 — 25 s _1 ; 

in the logarithmic x-axis = 0-01 is set to correspond to 

= 0, i.e. the standard recombination. Note differences be- 
tween the panels on y-axes scales. 



logical gas dynamic simulation, although current resolution 
with AM ~ 10 6 Mq seems not to be sufficient. 

However, with the upcoming low frequency interferom- 
eters (LOFAR and SKA) such modeling seems to be ex- 
cessive and unnecessary. Indeed, the absorption line pro- 
files and the spectral features from separate minihaloes lo- 
cated at z = 10 are of 1 to 5 kHz in width - c l ose to the 
resolution limit 1 kHz. As mentioned IXu et al. I (|201 lh the 
sensitivity corresponding to such a resolution requires enor- 
mously bright background sources, GRB afterglows and/or 
QSOs, and broad-band observations can be a better alter- 
native. In broadband observations the suppression of opti- 
cal depth in 21 cm caused by decaying particles manifests 
as a factor of 2-4 decrease of absorption in the frequency 
range v < 140 MHz where contribution from low mass mini- 
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haloes at z 10 dominates. Indeed, the number of haloes 
in the low-mass range 10 5 — 1O 7 M0 at z = 10 scales as 
n(M,z) ~ M _1 (the Press-Schechter mass function). The 
probability of minihaloes to intersect a line of sight is pro- 
portional to n(M, z) x (ar v i r (M)) 2 , resulting in a decrease of 
the number of strong absorption lines with W„ bs Si 0.3 kHz 
by factor of at least 2.5 for = 0.3 and more than 4.5 

for = 1. Such a decrease inevitably cancels the average 
21 cm forest signal from low mass haloes expected in future 
broadband observations. 



6 CONCLUSIONS 

In this paper we have considered the influence of de- 
caying dark matter particles on the HI 21 cm absorp- 
tion features from low mass minihaloes with mass M = 
10 5 , 1O 6 M , 10 7 M Q virialized at z vir = 10. We used a ID 
self-consistent hydrodynamic approach to study their evo- 
lution, and followed through the absorption characteristics 
from the turnaround to the virialization of minihaloes. We 
have found that 

• due to an additional heating brought by decaying par- 
ticles thermal and dynamical evolution of minihaloes in the 
presence of decaying dark matter shows cosiderable differ- 
ence from that in the model without particle decays (i.e. 
with the standard recombination scenario); 

• this additional heating strongly suppresses optical 
depth in the 21 cm line: it practically "erases" the 21 cm ab- 
sorptions from minihaloes with M = 10 5 — 1O 6 M0 even at a 
relatively modest (£ £ 0.3£l) ionization rate from decaying 
DM particles; the horn-like dependence of the optical depth 
found for minihaloes wi th M ~ 10 6 — 10 7 Mq in the standard 
recombination scenario IXu et al. I (|201ll ) almost disappears 
even at a lower ionization rate £ £ 0.1£l; 

• the equivalent width of the 21 cm absorption line de- 
creases significantly while £ increases, and the number of 
strong absorption lines W° hs ^ 0.3 kHz at z = 10 drops by 
more than 2.5 to 4.5 times depending on f; such a decrease 
inevitably erases the averaged 21 cm forest signal from low 
mass haloes (M ~ 10 6 — 1O 7 M0, i.e. the frequency range 
v < 140 MHz) in future broad-band observations. 
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